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Abstract. Disordered systems exhibiting exponential localization are mapped to anisotropic spin chains 
with localization length being related to the anisotropy of the spin model. This relates localization phe- 
nomenon in fermions to the rotational symmetry breaking in the critical spin chains. One of the intriguing 
consequence is that the statement of Onsager universality in spin chains implies universality of the localized 
fermions where the fluctuations in localized wave functions are universal. We further show that the fluctu- 
ations about localized nonrelativistic fermions describe relativistic fermions. This provides a new approach 
to understand the absence of localization in disordered Dirac fermions. We investigate how disorder affects 
well known universality of the spin chains by examining the multifractal exponents. Finally, we examine 
the effects of correlations on the localization characteristics of relativistic fermions. 

PACS. PACS-key Harper equation, localization - PACS-key Dirac fermions 
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Two-dimensional Ising model is one of the few exam- 
ples of exactly solvable many body systems. |2| The model 
exhibits a phase transition at finite temperature character- 
ized by universal exponents defining a universality class, 
the Onsager universality which describes the phase tran- 
sitions for anisotropic XY models. Interestingly, the On- 
sager universality also describes a quantum phase tran- 
sition at zero temperature, driven by quantum fluctua- 
tions, of one-dimensional quantum anisotropic XY spin 
chains in a transverse field. B These quantum models be- 
long to a small family of integrable Hamiltonians that 
have attracted both theoreticians as well as experimen- 
talists. Heart of integrability of these many body quan- 
tum spin problem is a mapping between spin and fermions 
which relates interacting XY spin chain with 0(1) sym- 
metry to the fermion Hamiltonian that are quadratic in 
fermions. The spin-fermion correspondence has proven to 
be extremely important also for the case of disordered 
magnetic field as the quadratic fermion Hamiltonian can 
be numerically diagonalized with extreme precision.^] Re- 
cent studies have shown that the a large variety of disor- 
dered quantum chains with 0(1) spin symmetry are still 
described by the Onsager universality. 
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In this paper, we show a new type of spin-fermion 
mapping where disordered fermions exhibiting exponen- 
tial localization are shown to be related to anisotropic spin 
chain in a disordered magnetic field at the onset to long 
range order (LRO). This correspondence, valid exclusively 
for exponentially localized systems, provides new insight 
into various issues relevant to disordered fermion as well as 
spin models. In particular we exploit well established uni- 
versality hypothesis of spin systems to make important 
statements about fermion problems. Firstly, the root of 



recently observed universality in localized fermions |lO|.[llj] 
is traced to the Onsager universality of the spin systems. 
Another interesting result is a correspondence between the 
relativistic and the nonrelativistic fermions in the presence 
of disorder. It is shown that the relativistic fermions can 
be viewed as the fluctuations in the exponentially local- 
ized solutions of the nonrelativistic fermions. This pro- 
vides a new approach for understanding the absence of 
localization in disordered Dirac fermions which has been 
the subject of various recent studies. (l^] We argue that 
the long range magnetic correlations provide mechanism 
for derealization of relativistic fermions thus obtaining a 
deeper understanding of the absence of localization. In ad- 
dition to obtaining intuitively appealing picture of some 
surprising results of disordered fermions, we also obtain a 
generalization of universality statement of the critical ex- 
ponents for disordered spin models. Finally, we examine 
how the correlations affect the localization characteristics 
and show the possibility of derealization of relativistic 
fermions analogous to the corresponding nonrelativistic 
case. 

Although the setting we describe is quite general in 
the context of disordered systems, for concreteness we 
will consider quasiperiodic disorder where the lattice prob- 
lem for nonrelativistic fermion is described by the Harper 
equation, Q 

lp n -l + 1p n +l + 2XV n 1p n = Elp n . (1) 

Here V n — cos(9 n ) where 9 n = (2na g n + <fi). The a g is 
an irrational number describing competing length scales 
in the problem and is a constant phase factor. Harper 
equation, in one-band approximation describes Bloch elec- 
trons in a magnetic field. This problem frequently arises 
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in many different physical contexts, every time emerging 
with a new face to describe another physical application. 
The problem is solvable by Bethe-ansatzQ: there is an 
algebraic Bethe-ansatz equation for the spectrum. It was 
shown that at some special points in the spectrum, e.g 
at the mid-band points, the Hamiltonian in certain gauge 
can be written as a linear combination of generators of 
a quantum group called U q {sl2)- It also describes some 
properties of the integer quantum Hall effect: it showed 
that Kubo-Greenwood formula for the conductance of any 
filled isolated band is an expression for a topological in- 
variant, and is an integer multiple of e 2 /h. It was further 
shown by Avron et al@ that it is a topological invariant 
that defines the first Chern class of the mapping of the 
Brillouin zone( a two-dimensional torus) onto a complex 
projective space of the wave functions. 

In contrast to usual Anderson problem describing lo- 
calized particle in a random potential, the Harper equa- 
tion exhibits localization-delocalization transition provided 
<j g is an irrational number with good Diophantine prop- 
erties (ie. badly approximated by rational numbers). This 
transition at A = 1 has been characterized with singular 
continuous states and spectra. Richness and complexity 
of the critical point describing localization transition has 
been studied in great detail by various renormalization 
group approaches. [^[|^] Recently, it was shown that mul- 
tifractal characteristics continue to exist beyond the criti- 
cal point throughout the localized phase. This hidden 
complexity of the localized phase is brought to light after 
one factors out the exponentially decaying envelope. The 
localized wave functions with inverse localization length 
= log(X) is rewritten as, [fill 
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where 7 = £ . The tight binding model(tbm) describing 
the fluctuations rj n in the exponentially decaying envelope 
is given by the following pair of equations, 



r l r n+l + ^rf n _ l +2XV n r l r n = Erf n 
2XV n rf n = Eri 
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Fig. 1. Absolute value of the fluctuations for Harper equation 
for E = states with <j> = .25 and a g = (-\/(5) - l)/2. At the 
Fibonacci sites, we see period-6 behavior(period-3 in absolute 
value): \r\w m \ = \rjF m+3 \ which is independent of A. 



a transverse magnetic field, given by the following spin 
Hamiltonian. 

H = + e 7 «+i + 2AK<]. (4) 

The Cfc, k = x,y are the Pauli matrices. The e -7 and 
e 7 respectively describe the exchange interactions along 
the x and the y directions in spin space. Using Jordan- 
Wigner transformation, the interacting spin problem can 
be mapped to non-interacting spinless fermion problem[|| 
where fermions are the quasiparticle excitations of the spin 
chain obeying the following coupled equations, 



Here rf and rf respectively describe the fluctuations to 
the right and to the left of the localization center. An 
exact renormalization scheme JlO| showed that these fluc- 
tuations exhibit universal features(see Fig. 1) described by 
the strong coupling fixed point. Hence, localized phase is 
characterized by universal fractal characteristics described 
by A — > 00 limit of the equation. These results were further 
confirmed by rigorous mathematical analysis. [jllj Hence 
the strength of the disorder which determines the localiza- 
tion length can be factored out making universal aspects 
transparent. 

A new revelation that provides intuitive understand- 
ing of the strong coupling fixed point of Harper is ob- 
tained by relating the fluctuations described by (3) to the 
anisotropic spin chain at the onset to LRO. It turns out 
that the equations (3) for E = describe quasiparticle 
excitations of a critical anisotropic XY spin- 1/2 chain in 
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2\V l r ] } = Er# 
4 2AV^f = Er,] 



(5) 



At the onset to LRO, the excitation spectrum becomes 
gap-less and hence for the E = 0, the massless mode, 
the above two equations are degenerate coinciding with 
the equation (3). Therefore, the massless excitations of 
the critical anisotropic spin chain describe the fluctuations 
in the exponentially localized excitations of the isotropic 
chain where the anisotropy parameter is related to the 
localization length. 

Therefore, the localization length like the anisotropy 
parameter is an irrelevant variable and hence the state- 
ment of strong coupling universality of the Harper equa- 
tion is synonymous with the statement that anisotropic 
spin chain is in the universality class of the Onsager uni- 
versality. In short, we establish an equivalence between 
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Ising fixed point of the anisotropic spin chain and the 
strong coupling fixed point of the Harper equation and 
thus provide a simple interpretation of the strong cou- 
pling universality of the Harper equation. Furthermore, 
by relating anisotropy to the localization length, we ob- 
tain a new picture of the localized phase: the role of the 
strength of the disorder is similar to that of the strength of 
the anisotropy, which is an irrelevant parameter for renor- 
malization flow. It should be noted that the spin-fcrmion 
mapping provides a new method to determine the local- 
ization length of the tbms in the presence of disorder. The 
expression for the localization length, 7 = log(X) can be 
viewed as the relation describing the critical point of the 
spin chain in a disordered field. 

We next show that the fluctuation rj n obey Dirac equa- 
tion for zero energy states. In the long wave length limit, 
the equations(3) reduce to the Dirac equation. We replace 
n by x and write r\ n ±\ — e ±lp r)(x), where p is the momen- 
tum canonically conjugate to x. The equation (4) for the 
fluctuations for E = state can be described by the fol- 
lowing non-Hcrmitian Hamiltonian Hfi uc and its adjoint, 
Hfiuc = e~^e lp + eie~ ip + 2XV(x). In the limit {p -> 0), 
the system for E — reduces to the Dirac equation, 



[ga x p - i(2XV{x) + 2cosh(j))a y ]r)(x) = (6) 

where ij(x) is a two-dimensional spinor r){x) = (rj l (x), T] r (x)) 
It is interesting that the two-component structure of Dirac 
spinor arises naturally when we consider fluctuations about 
exponentially localized wave functions. Here g = 2sinh(-y) 
is the velocity of the Dirac fermions while the mass of the 
Dirac fermions m(x) = 2((XV(x) + cosh{^f)). Therefore, 
on a lattice, the Dirac fermions with disordered mass are 
the fluctuations of the nonrelativistic localized fermions. 
This would imply the absence of exponential localization 
for relativistic fermions. The defiance of localization by 
relativistic fermions has been the subject of various stud- 
ies and our analysis provides a simple way to understand 
this intriguing result. 

Next, we address the question whether the strong cou- 
pling fixed point which describes localized phase of Harper 
equation, the critical Ising model and the Dirac fermions 
with quasiperiodic disorder, implies universal multifractal 
exponents. We compute the f(a) curve( Fig. 2) describ- 
ing the multifractal spectrum associated with the self- 
similar wave function or the inverse participation ratios 

P, P(q,N) = ~ N-^\ a=% and f{a) = 

aq — r(g) The free energy function r(q) and its Legen- 
dre transform f(a) were found to be A independent only 
for for positive values of q and hence only left half of the 
f(a) curve is universal. Therefore, for quasiperiodic spin 
chains at the onset to LRO, scaling exponents for only 
positive moments of the participation ratio are universal. 
This can be viewed as a generalization of the universality 
statements of periodic spin chains to disordered spins. 
Finally we investigate the role on correlations on the 




Fig. 2. Numerically obtained f(a) curves for A — > 00 ( solid 
curve) and A = 1.5 (lines with crosses). 



For quasiperiodic disorder, dimer-type correlations can be 
introduced by replacing 6 n = 2na g n in V(6 n ) by the iter- 
ates of the supercritical standard map, describing Hamil- 
tonian systems with two degrees of freedom. [jl5| 

n +i + n -i - 29 n = - — sin(27T0 n ). (7) 

We use iterates that describe golden-mean cantorus(the 
remanent of the KAM torus beyond the onset to global 
stochasticity) which has been shown to exhibit dimer-type 
correlations, and leads to Bloch-type states for the non- 
relativistic fermions. Here, we will confine to the Ising 
limit( can be obtained from (5) using 7 — > 00 and rescaling 
the parameters), described by, 



Vn+l 

vl-i 



2Xcos(27T0 n )n 1 n = Erft 
2Xcos(2n0 n )r 1 2 n = En\ 



(8) 



We determine the critical A ( the threshold for the onset 
to magnetic transition) as a function of K, the nonlinear- 
ity parameter of the two-dimensional map.Q The local- 
ization characteristics of the massless mode of the Ising 
model are studied using an exact RG methodology M. 
In this approach, quasiperiodic models such as eqn (7) 
with golden-mean incommensurability are decimated to 
a renormalized model defined only at the Fibonacci sites. 
Rcnormalization flow describing renormalized couplings at 
the Fibonacci sites provides an extremely accurate tool to 
localization characteristics of massless spin excitations which distinguish extended, localized and critical states. Trivial 
obey Dirac equation. The fact that correlations would re- fixed points of the RG describe extended states while crit- 
sult in derealization as originally shown by a random- ical states correspond to nontrivial fixed points. As shown 
dimer model (lj] is an important result in localization theory. [0Fig. 3, nontrivial 6-cycle ( which also corresponds to six- 
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Fig. 3. (a)Critical A as a function of A" for the Ising model, (b) 
RG 6-cycle showing the variation in the renormalized coupling 
at Fibonacci iteration of the RG flow for E = 0. 

cycle of the wave function rj n as shown in Fig. 1) degener- 
ates to trivial fixed points at certain special parameter val- 
ues. The origin of these trivial fixed points, has been traced 
to a hidden dimer in the quasiperiodic iterates describing 
the golden-cantorus . p5| At these points, the relativistic 
massless mode of the Ising model is ballistic. Therefore, 
relativistic fermions may become delocalized analogous to 
the nonrelativistic case due to dimer- type correlations. 

The Fig. 3 shows an interesting interplay between the 
magnetic transition and the ballistic transition due to 
dimer-type correlations: the ballistic transitions where the 
relativistic mode is propagating seems to be sandwiched 
between two peaks corresponding to strong enhancement 
of ( possibility divergent ) strength of the inhomogeneous 
field needed for the onset to LRO. This phenomena again 
confirms the view that spin-fermion relationship may be 
an extremely useful means to understand the richness and 
complexity underlying a variety of new phenomena in dis- 
ordered systems. 

One-dimensional quantum spin chain in a transverse 
field at the onset to LRO describes two-dimensional lay- 
ered Ising model. Therefore, our study relates universal as- 
pects, described by Onsager universality, of two-dimensional 
Ising model to the universal aspects of the two-dimensional 
Bloch electron problem described by the Harper equation. 
This paper establishes a relationship between two impor- 
tant systems where geometry and integrability are of cen- 
tral importance. We believe that our results are valid for 
a variety of disorders including a large class of pseudoran- 
dom as well as random cases. 

Finally, it should be noted that spin-fermion corre- 
spondence is also valid for time-dependent models. Kicked 



Harper model that has been extensively studied in quan- 
tum chaos literature also describes XY spin chain 
in a periodically kicked inhomogeneous magnetic field. Jl^| 
By exploiting the spin-fermion correspondence, various 
new results in the kicked Harper model can provide a new 
way to understand many surprising results. [|l9| One of the 
interesting results is the independence of the critical ex- 
ponents with respect to the discommensuration parame- 
ter a g in the limit u g — > 0. This defines a new aspect of 
the usual universality statements for the spin systems at 
the onset to magnetic transition and hence broadens the 
concept of universality to include disorder as well as time 
dependence. 
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disordered Dirac fermions. 



References 

1. Web address: physics.gmu.edu isatija 

2. L. Onsager, Phys Rev 65, 117 (1944). Also see H. E. Stan- 
ley, " Introduction to Phase Transitions and Critical Phe- 
nomenon", McGraw-Hill, 1978. 

3. E. Lieb, T. Schultz, and D. Mattis, Ann Phys. (N.Y) 16, 
407 (1961). 

4. J. M. Luck, J . Stat. Phys 72 (1993) 417; Eurpphys Lett. 24 

(1993) 359. 

5. LI. Satija, Phys Rev B49, 3391 (1994). 

6. P.B. Wiegmann and A.V. Zabrodin, Phys Rev Lett 72, 1890 

(1994) . 

7. J.E.Avron, R. Seiler and B. Simon, Phys. Rev. Lett. ,51, 
1983. 

8. S. Ostlund and R. Pandit, Phys. Rev. B 29, 1394 (1984); S. 
Ostlund, R. Pandit, D. Rand, H. J. Schellnhuber, and E. D. 
Siggia, Phys. Rev. Lett. 50, 1873 (1983). 

9. J. A. Ketoja and I. I. Satija, Phys. Lett. A 194, 64 (1994), 
Physica A, 219, 212, (1995). 

10. J. A. Ketoja and I. I. Satija, Phys. Rev. Lett. 75, 2762 

(1995) ; 

11. B.Mestel,A.Osbaldestin and B. Winn;"Golden mean 
renormalization for the Harper equatiomthe strong coupling 
fixed point" , preprint. 

12. L. Balents and M. Fisher, Phys Rev B, 56, 12970 (1997). 
Castillo, C. Chamon, E. Fradkin, P. Goldbart and C. Mudry, 
Phys Rev B, 56, (1997) 10668; 

13. P. G. Harper, Proc. Phys. Soc. London A 68, 874 (1955); 
B. Simon, Adv. Appl. Math. 3, 463 (1982); 

14. H-L. Wu, W. Golf and P. W. Phillips, Phys Rev B 45 
(1992) 1623. 

15. J. Ketoja and I. Satija, Phys Rev B, 59, 9174 (1999); 

16. P. Leboeuf, J. Kurchan, M. Feingold and D. P. Arovas, 
Phys.Rev.Lett 65, 3076 (1990). 

17. T. Prosen, I. I. Satija and N. Shah, Phys. Rev. Lett., 87, 
066601-1 (2001) 

18. T Prosen and I Satija (unpublished). 

19. "h — * in kicked Harper": Reassurances and Surprises" , 
I Satija and T Prosen ( preprint). 



